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On the Solutions of Linear Homogeneous Difference 

Equations.* 

By E. D. Cakmichael. 



Introduction. 

Let us consider the linear homogeneous difference equation of order n, 

F(x+n) + a 1 (x)F(x+n—l) +a 2 (x)F(x+n—2) + + a n (x)F(x) = 0, (1) 

in which the coefficients are analytic at infinity or have poles there. By means 
of a transformation of the form 

we arrive at the new equation, 

f(x+n) +a x (x)f{x+n-l) + a 2 (x)f(x+n-2) + .... +a n (x)f(x) =0, (2) 

where 

a k (x)=(x+n— k)^ x + n - k >(x+n)-^ x+n) a k (x), k = l, 2, , n. 

It is obvious that integers (i exist such that each of the functions a k (x) is an- 
alytic at infinity. For the value of [i we choose the least integer such that the 
functions a k (x) are all analytic at infinity.! Let us write 

a k (x) = a m + -2SL + f". + . . . ., k = l, 2, . . . ., n, \x\>R. (3) 

We shall confine our attention to the case in which the constants a 1Q , a 20 , 
. . . . , a n0 are such that the characteristic algebraic equation \ 

a"+a 10 <x M - 1 + a 20 a' , - 2 + .... +a n0 =0, (4) 

has its roots a x , a 2 ,...., a n , different from each other and from zero. In this 

* Read before the American Mathematical Society at Chicago, December 29, 1914. 

t It is easy to see that x—^ak (x) is analytic at infinity for every value of k from 1 to n and that 
/» is the least integer for which this is true. We shall say that the value of aj-0*ok (x) at infinity is Uka . 
Comparing with equation (3) we see that a k d=e-f'''ako- 

$ It will be said that this is the characteristic algebraic equation associated either with (2) or 
with (1). 



186 Cabmichael : On the Solutions of Linear 

case, as is well-known, equation (2) has n formal power series solutions of the 
form * 

/,(*) = a?0*(l+ -^ + ^ +....), »=1, 2, . . . ., n. (5) 

These are in general divergent. 

In this paper I shall show how the formal solutions (5) afford a means by 
which equation (2) may be separated into two members so that the method of 
successive approximations, when applied in its usual form, will lead to an actual 
solution of (2) and by aid of this single solution shall develop means by which 
the most important fundamental systems of solutions of (2), and hence of (1), 
are found. 

This paper contains no new results, so that its novelty lies entirely in the 
new methods employed. The theorems demonstrated are essentially equivalent 
to those obtained by Birkhoff in the Transactions of the American Mathemat- 
ical Society, Vol. XII (1911), pp. 242-284, and hence are more complete than 
the similar ones previously demonstrated by me and published in the volume 
just referred to, pp. 99-134. The method of the present paper is more closely 
related to that of my previous paper than to that of Birkhoff's later memoir. 
It shares with both of them the desirable property of being a direct means for 
attaining the end in view ; but it is essentially different from either of them and 
in several respects is to be preferred to them. I have found it convenient to 
use the Birkhoff modification of the integral earlier employed by me, since this 
modified integral serves the purpose of securing a rather simpler fundamental 
system of solutions, especially as regards the asymptotic character of the func- 
tions in the solution. 

For other methods of dealing with the problem of difference equations see 
Norlund, Acta Mathematica, Vol. XL, pp. 191-249, and the papers to which he 
refers. 

§1. Formal Solutions by Successive Approximation. 
Let us denote by l{(oo) the functions 

k{x) =a?»*(l+ C f + % + .... + £•), i=l, 2, . . . ., n. 

Form the equation 



* The value of the constants m , ca , on , . . . . may be determined by substituting the value of ft (w) 
from (5) into (2) and directly reckoning out the constants so that the resulting equation shall be a formal 
identity in x. 
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g x {x+n) g^x+n—1) g t (x) 

Ijix+n) l 1 (x-\-n—l)....l 1 (x) 
l % (x-\-n) l 2 (x-\-n — 1) . . . . l 2 {x) 



:0. 



(6) 



l n {x+n) l n (x + n—l) l n (x) 

This is obviously a linear homogeneous difference equation of order n at most ; 
and it is of order n if the minors of g x {x-\-n) and g x {x) are neither identically 
zero. These minors are D(x) and D(x+1) respectively, where 

l^x+n—1) l^x+n—2) l^x) 

l 2 (x-\-n — 1) l 2 (x-\-n — 2) . . . .l 2 (x) 



D(x) 



(7) 



l n (x+n— 1) l n (x+n— 2) l n (x) 

Divide the first row of this determinant by o^af 1 , the second row by o^x^, . . . ., 
the last by a£a/"\ Each element of the resulting determinant is analytic at 
infinity; and hence the function D(x) represented by this resulting determinant 
is analytic at infinity. It is easy to see that the value of D{x) at infinity is 



Z)(oo) = 



a? -1 
aT 1 



ar 2 - 
ar 2 . 



OCl 

a 2 



• a„ 



or 



D(ao) = H(a r — a s ), r<s. 



Since the numbers a x , a 2 ,...., a n are all different by hypothesis, we see that 
D(x) is not identically zero; in fact, it is different from zero at infinity. 
Hence, D(x) is not identically zero: and therefore equation (6) is of order n. 
In the determinant in equation (6) divide the second row by afa" 1 , the 
third row by o^x^ 2 , . . . . , the last row by a*^", and expand the resulting deter- 
minant in terms of the elements of the first row. The coefficients of the quan- 
tities g^x-j-n), g x (x-\-n — 1), . . . ., g t (x) are obviously all analytic at infinity, 
that of g-tix+n), namely, D(x), being different from zero at infinity. Hence 
we have from (6) an equation of the form 

g^x+n) +\(x)g 1 {x+n-l) +^(x) gi (x+n-2) + .... +?, n (x)g 1 (x) = 0, (8) 

in which \{x), Xz(x), . . . ., l n (x) are analytic at infinity. This equation has 
the solutions l x {x), k(x), ....,l n (x), as we see readily from equation (6). 
Moreover, these form a fundamental system of solutions, since D(x) is not 
identically zero. 
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By means of the fact that equation (8) has the solution lt{x), h(x), . . . ., 
l n (x) and the relation which exists between l^x), i = l, 2, . . . ., n, and the for- 
mal solution (5) of equation (2), it may be shown by direct computation that 



X 



X" 



X 



that is, that the expansion of ^(x) in a descending power series in x coincides 
with that of a t (x) as far as the term in l/x m+1 . 
Now, write equation (2) in the form 



2 \{x)f{x+n-k) = 2V k (x)f(x+n-k), 7*{x)=l. 
Then the function ty k (x) may be expanded as follows: 

k = l, 2, ,n. 



jjf / „\ Pk, ro+2 i Pk,m+3 , 



Form the system of equations 

2 A, 4 (fl5)0i(a5+w— &) =0, 

fc=0 
n re 

2 a fc (a0# 2 (« + w— A;) = 2 W k (x)g 1 (x + n—k), 



(9) 



(10) 



&=0 



ft=l 



2 X t (as)^(aj+n— &) = 2 ^ k (x)g i _ 1 (^ + n— k), 

fc=0 fc=l 



(11.) 



By means of these, used as equations for successive approximation, we shall 
be able to obtain formal solutions of equation (9) or (2). 

The first equation in (11) is homogeneous and has the fundamental system 
of solutions h{x), l 2 (x), . . . ., l n (x), as we saw above. The remaining equa- 
tions in (11) are of the form 



2 X k (x)g(x+n— k) =*i(x). 



fc=0 



(12) 



We shall next find two formal solutions of this typical equation. Our method 
is the classic one of variation of parameters. We assume a solution of the form 



g(x)= 2t k {x)l k (x), 



k=l 



(13) 



and determine the functions ^(x), t 2 (x), . . . ., t n (x) in a convenient way so 
that they shall formally satisfy equation (12). Thus we may write 
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g(x+l) = 2t k (x)l k (x+l), if 2 At k (x)l k (x + 1) = 0; 



fc=i 



g(x+2)=2t k (x)l k (x+2), if 2^(a)l,(*+2)=0i 



g(x+n—l)=2t h (x)l k (x+n-l), if 2 A*»(aj)^(»+n-l) =0; 



j;=i 



fc=i 



g(x+n) — 2 <*(»)^(af+n) + 2 At k (x)l k (x+n). 
k=i k=i 

So far we have imposed n — 1 conditions on the n functions h(x), t 2 (x), . . . ., 

t n (x). We obtain an «-th condition, namely, 

n 

2 At k (x)l k (x+n) =7i (x), 

k=l 

by substituting the preceding values of g{x), g(x-\-l), . . . ., g(x+n) in equa- 
tion (12), remembering that l x {x) ,l 2 {x) , . . . . ,l n {x) satisfy the equation obtained 
from (12) by replacing ri(x) by zero. 

The n equations of condition now placed upon the functions t{x) may 
readily be solved by determinants for the quantities At k (x), k = l, 2, . . . ., n. 
Thus we have 

At k {x) = (-iy +k v (x)M k (x + l)/M(x + l), k = l, 2, 

where 

hi*) h(x) l n (x) 

h(x+l) h(x + l) ....l n (x + l) 



.,n, 



:w) 



M(x) = 



^(x+n— l)l 2 (x+ n— 1) .... l n (x + n— 1) 
and M k (x) is the minor of the element l k (x+n — 1) in M(x). 
Now, the function 

M(x)/l 1 (x)l 2 (x) l n (x) 

is analytic at infinity and is different from zero at that point, as one may see 
readily by considering its value in the form of a determinant obtained from 
that for M(x) by dividing the columns in order by l x (x), l 2 (x), ... ,,l n (x); 
for each element of the resulting determinant is analytic at infinity and the 
determinant has at infinity the same value as the determinant D(x) consid- 
ered above. Likewise it may be shown that the function 
M k {x)/l x {x) h-i{x)l k+1 (x) l n (x) 

is analytic at infinity and is different from zero at that point. Hence, we may 

write 

(~iy+*M k (x)/M(x) =A k (x)/l k (x), 

24 
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where A k (x) is analytic at infinity and is different from zero at that point. 
Hence, we have from (14) 

At k (x)=n(x)A k (x+l)/l k (x+l), k = l, 2, . . . ., n. (15) 

Equations (15) have the following two formal solutions : 

CO 

t k {x)— — 2,yi(x+i)A k (x+i + l)/l k (x+i + l), k=l, 2, , n; 

t h {x) = X,ri(x — i)A h {x — i + l)/l k (x — i + 1), k = l, 2, . . . ., n. 

Substituting these values in (13), we have the following two particular 
formal solutions of equation (12) : 

g + (x)= — 2 i V (x+i)A k (x+i+l)l k (x)/l k (x+i+l); (16) 

co n 

g-(x) = I, 2ri(x—i)A k (x—i+l)l k (x)/l k (%—i+l). (17) 

These formal solutions will be denoted by 

S,(>l) and T x {v,), 
respectively. 

Let Px(x), p 2 (x), ... ., p n (x) be n arbitrary periodic functions of x of 
period 1. Then for the formal solutions of the successive equations in (11) 
we may conveniently select either of the following sets : 

n 

gt{x)= 2p k (x)l k (x), 

k=l 
n n n 

g£(x) = Xp k (x)l k (x)+S x \ 2%(x) 2 p k (x)l k (x+n— i t )\, 

k-l 4=1 k=l 

n n n 

gt{x)= 2p k (x)l k (x)+S x \ 2 %(x) X p k (x)l k (x+n—i 1 ) \ 

k=l ii=\ k=l 

+ S x [2%(x)S x+n _ k \iv il (x)Zp k (x)l k (x+n~i 1 )\], 

i 2 =l ii — 1 &=1 



gi(x) = 2p k (x)l k (x), 

fc=l 
n n n 

g;(x) = 2 p k (x)l k (x) + TJ 2 %(x) 2 p^x^x+n-^) j , 

fc=i »,=i fc=i 



g;(x) = 2 Pk(x)h(x) +T X \ 2 % t (x) 2 p^x^ix+n-i,) 



n n 



+ T x [2%,(x)T x+n _ i J X%(x)Sp k (x)l k (x+n-i 1 )\], 
i g — i £1=1 fc=i 
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We are thus led to the following two forma] expansions : 

f + (x)=2p k (x)l k (x)+S x \ i%(x)ip k (x)l k (x + n-i 1 )\ 

+8 x [i%(x)8 x+n _ k \ 2 %(x)ip k (x)l k (x+n-i 1 )\] + ...., (18) 

t(x)=£ Pk (x)l k (x)+T x \ i% t {x)ip k {x)l k (x+n-H)\ 

fc=l i 1= l fc=l 

n n n 

+ T x [Z%(x)T x+n _ h \ZV h (x)2p k (x)l k (x+n-i 1 )\] + (19) 

It is easy to verify that these functions afford formal solutions of equation 
(9). It is sufficient to this purpose to make a direct substitution of the expan- 
sions into the equation, reducing the result by means of the relations 

n n 

X \{x)S x+n _ k {Yl) = V}, 2 \{x)T x+n _ k (Yl) — Y]. 
ft=0 k=d 

Each of these solutions is in general illusory on account of the divergence 
of certain series contained in the terms of the expansions. Let <x r and a, be 
any two a's such that |a,|>|a s |. Put p { (x)=0, i^r; p r (x)=l. Then the 
second term in the series in (18) takes the form 

S x \ i%(x)l r (x+n-i l )\=-2 i%(x+i)A k (x+i+l) 

i'i=l «'=0 k=l 

l r (x+n+i—i 1 )l k (x)/l k (x+i + l). 

It is easy to see that the series in the second member of this equation is not in 
general convergent. A similar discussion may be made of the second term in 
the series in (21) ; it turns out that the series denoted by this term is not in 
general convergent. 

§ 2. Two Actual Solutions. 

From the last result of the preceding section we see that the series in (18) 
and (19) do not in general represent actual solutions of equation (9). We 
shall now show that an appropriate choice of the periodic functions Pi(x), 
p 2 ( x ),...., p n (x) will lead to actual solutions of this equation. 

First, we shall consider the case of the series in (18). We assume that 
the notation for the roots a x , <x 2 ,...., a„ of the characteristic algebraic equa- 
tion is chosen so that 

I«i|<k2|<....<|a„|. (20) 

Then put 

P X {X) -1, p 2 (x) —0=p s (x) = =Pn(%)- 
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Denote the resulting formal solution (18) by ft{x). Then we may write 

ft(x)=l,(x)+8 x \ £«P<»*iOH-»-*.){ 

+S x [£%(x)S x+n _ i J ZV^kfr+n-H) }] 

+ #7 £ *«.(*) A.+.-J £ ^»<Wd 2 *.•»*! (&+«-%) | J 

\i 3 =l i 2 =l ii=l 
+ 



(21) 



It is convenient to introduce the term right D-region to denote a region of 
the complex plane defined as follows : Consider .the curve made up of a semi- 
circle C with center at the point zero and terminated at each end on the axis of 
imaginaries, its position being to the right of that axis and its extremities being 
A and B, together with the straight lines Acq and Deo drawn from the points 
A and B, respectively, to infinity in a direction parallel to the axis of reals and 
lying in the left half -plane. This curve <x> ACBoo divides the plane into two 
parts, one containing the point and the other not containing this point. The 
part not containing the point is called a right D-region. It is understood 
that this D-region does not contain the point infinity. 

Now it is clear that there exists a right D-region, call it D x , and constants 
M 1} M 2 , M s such that each of the following inequalities is true so long as x is 
inD i: 

<M 1 |*|- m - 2 , 



A k (x+i + l)l k (x) 

l k (x + i + l) 
l x {x-\-n — k) | 
h(x+i) | 



<M 2 \a k - i \, 

<tf,|*i(aO|, 

<2\h(x)\-\ 



ai| 



k = l, 2, 



n. 



(22) 



It will also be assumed that D x is chosen so that the functions ^ h {x), A k (x), 
\l k (x) \~ % are analytic in D x . 

The first term in the second member of equation (21) is a function which 
is analytic throughout the finite plane except at x=0. The second term denotes 
a series. It is obvious that 



\8 X \ ZV^hix+n-iJWSSsl X%(x)l 1 (x+n-i 1 )\, 

where S x denotes the series obtained from 8 X by replacing each term in the 
latter by its absolute value. If we confine x to the region D x , and employ 
inequalities (22) and (20), we come readily to the following relations among 
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infinite series, these relations being valid term by term : 

S,\ i%(x)l 1 (x+n-i 1 )\\<2M 3 \l 1 (x+i)\.nM 1 \x+i\- m -*-M 2 i |ajp| 

< i | h {x) | • 2n 2 M AM 3 ! oo+i \ ~ m ~\ ( 23) 

If we further confine x to lie in some closed region 8 in D 1 , then it is clear 
that a constant M exists such that the last series above is term by term less 
than the series of constant terms 

XM(i+l)- m ~ 2 . 

This series being convergent, it now follows by the theorem of Weierstrass that 
the series denoted by 

S,\ i^Whix+n-iJl (24) 

is uniformly convergent in 8. Furthermore, each term of this series is analytic 
in S. Hence, the sum of the series is a function of x which is analytic through- 
out S. It is therefore analytic at every point in B 1 , since S is any closed region 
in B x . 

By means of (23) we shall now obtain more convenient inequalities gov- 
erning the character of the function in (24). "We have 

\8.\ *k% x {x)l x (x+n-i x ) | I <2n 2 M 1 M 2 M 3 \l 1 (x) |2 \x+i\~ m ~ \ 

If x=u + vV — 1 and u and v are real it is easy (cf. Birkhoff, loc. cit., p. 248) 
to show that we have 



=o |a5+*|" + " \x\ 

l 1 2n I I ^ -, 



<=o \x+i\ m+2 \v\ m+1 

Hence, if we assume that the semicircle used in defining D 1 has a radius greater 
than unity, we have the following inequalities valid in D x : 

IS S21P lx)l(x+n-i)\\<\ 2nn * MMiM s\ h{x) \'\ X \~ m ^ U >°' 1(25) 

It is obvious that the regions in which the two inequalities in (25) are sepa- 
rately valid overlap and that the two regions make up the whole of D 1 . 

We may proceed in like manner to the study of the function defined by the 
third term in equation (21). It is obvious that the series denoted by S x in the 
expression 
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S x [Z% 2 (x)S x+n _ i2 \ SVtWkfr+n-iJtt (26) 

i s =l 4=1 

is term by term equal to or less in absolute value than the series denoted by 
S x in 

S x [£%(x)S x+n _ i2 \ £«MaOIi(*+»-h) }]• (27) 

i 2 =l »'i=l 

It is now convenient to make a separation into two cases. 

Let us consider first the case in which «>0. By employing (25) we may 
readily see that the series denoted by S x in (27) is term by term less than the 
series 

S.[ £%(x)2nn 2 M 1 M 2 M 3 \l 1 (x+n-i 2 ) | • \x+n-i 2 \~ m ~\ 

and that this in turn is term by term less than the series 

S.[2wn l M 1 M^l^M 1 \x\- m ^M t \l 1 (x) | • |a?|— 1 ]. (28) 

From this point we may proceed as in the discussion of the second term 
of the series in (21). Thus we may show that the series denoted by S x in (26) 
is uniformly convergent in any closed region 8 lying in D x and having no 
point to the left of the axis of imaginaries, and therefore that (26) denotes 
a function which is analytic at every point in D a and not to the left of the axis 
of imaginaries. Moreover, through an examination of (28) it is easy to see 
that 

\S,[ £ %(x)S x+n _ h \ £ % (aO?t(»+»-h) [ J I 

<(2nn 2 M 1 M 2 M 3 y\l 1 (x)\-\x\- 2m - 2 , u>0. (29) 

For the next discussion we confine x to D 1 and to that part of T> x in which 
\v\ >1. Then the second inequality in (27) becomes available. Making use 
of this and proceeding by the method just employed in the preceding paragraph, 
we may show that the function (26) is analytic at every point in the region in 
consideration and that moreover the following inequality is satisfied : 

\S x [£%(x)S x+n _j£v il (x)l 1 (x+n-i 1 )\]\ 

< (4:7tn i M 1 M 2 M 3 ) 2 \ h(x) | • | v | - 2 "- 2 , | v | >1. (30) 

It is now easy to see that the method used above may be further employed 
in an investigation of the successive terms of the series in (21) and that the 
following results will thus emerge: Every term of the second member in (21) 
represents a function which is analytic at every point in D ± ; this series is term 
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by term less in absolute value than each of the following series : 



\hw\+\hw\^+\h(*)\(^y+ 



\h(x)\ + \l 1 (x) 



2M 



\m+l 



+ 1*1 0*0 I 



2M 



\m+l 



+ 



M = 2nn 2 M 1 M^U, 
, \v\>l, 



u >0 ; 



(31) 



| v | "*" ri ' ■ \ | v | 

each inequality being valid for that part of T> x for which the corresponding 
adjoined condition «>0 or \v\ >1 is satisfied. 

Let us put on the region D x a further restriction consistent with those which 
we have previously placed upon it, namely, that it shall have the property that 
for any x in D lf where x=u+vV—l, the following two inequalities shall be 
true: 

M\x\- m - l <1 when m>0; 2M|v|- m - 1 <l when |t>|>/8, 
/3 being an appropriately chosen constant. Now let 8 be any closed region 
lying in the modified D t . Then, from the result associated with (31) above, 
we see readily that the series in (21) is uniformly convergent in S ; whence it 
follows that the sum fi(x) is analytic at every point in the region D lf since the 
terms of that series are analytic at all such points. Moreover, it is clear that 
ft(x) satisfies the following inequalities, valid in D x : 

M 



\tt(x)-h(x)\<\h(*)\ 

\ft{x)-l 1 {x)\<\l 1 {x)\ 



X 



m+i 



If,' 



M>0 ; 



2M 



I •»+!. 



>fc 



(31') 



| v | "' Ti — 2M 

as one sees readily by transposing the first term in the series in (21) to the 
first member of the equation and comparing the resulting second member term 
by term with the series in (31), exclusive of the first term in each, and finally 
summing the last series (with first term removed) . 

From the last inequality we have readily the following limits : 



lim x m 



lim v m 



ft(x) 
l x {x) 



— 1 



= 0, 
=0, 



(32) 



the first being valid for any approach of x to infinity in the right half-plane, 
and the second being valid for the approach of the real variable v to either 
+ co or — oo . Since v= |a/|sin(arg#), it is easy to see from these two limits 
that 

7i + (a0 



lim x m \ 



f/i + (*0 A. 

\h(x) J 



:0, 



(33) 
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provided that x approaches infinity in any way so as to remain always without 
any sector determined by two rays from zero to infinity and including between 
them the negative part of the axis of reals. 

By actual substitution it may now readily be shown that the function fi(x) 
satisfies equation (2) for all x in I) 1 . By means of the equation 

f(x)= * f( x+n )-°^f{x+n-l)-....-^=^f(x+l), 

a n (x)' y a n (x)' a n (x) 

which is obtained from (2) in an obvious manner, one may further define the 
solution ft(x) in the part of the finite plane exterior to D x . It may thus be 
seen that the extended solution f?(x) is analytic throughout the finite plane 
except at the singularities of the functions 

1 Ms) g«-i(aQ (34) 

a„(x)' a n (x)' ' a n (x) 

and points congruent to them on the left. (A point a — i is said to be congruent 
to a on the left, if i is a positive integer; similarly, we say that a+i is con- 
gruent to a on the right, if i is a positive integer.) It is further obvious that 
the singularities of ft{x) are poles in case the functions (34) are rational and 
that the complete set of numbers each of which is the order of a pole of /f (x) 
is a bounded set. 

In obtaining the solution whose properties we have just developed, we 
began with the series (21) which at first was only known to have formal validity 
and by means of it arrived at the actual solution ft(x). Similarly, one might 
start from the formal solution 

f~(x)=l n (x)+T x \i%(x)l n (x+n-i 1 )\ 

+T X [ 2 %(x)T x+n _ ia \ 2 tp^WKix+n-i,) | ] + . . . ., (35) 

obtained from (19) in a way similar to that by which (21) was obtained from 
(18), and show that we are thus led to an actual solution f~(x) of (9). There 
would be no essential modification in the style of the argument. Instead of a 
right D-region we would employ a left D-region, such a region being defined 
by saying that its points are obtained from those of a right D-region by reflec- 
tion through the axis of imaginaries. In fact, practically the only modification 
of the argument necessary is that which comes of an interchange of the roles 
of the right and left sides of the plane. 

Combining the principal conclusions which we reached in the treatment of 
equation (21) and those which will emerge from a discussion of (35) in the 
manner indicated, we have the following result : 



Homogeneous Difference Equations. 197 

Preliminary Theorem. If the characteristic algebraic equation associated 
with the difference equation (2) has its roots a x , a 2 , . . . ., a re different from 
each other and from zero and these are ordered so that |a a | <|a 2 |<. . . . <|a„|, 
then this difference equation has two particular solutions ft{x) and fn(x), 
possessing the following properties : 

1) The solution ff(x) is analytic throughout the finite plane except at 
the singularities of the functions (34) and points congruent to them on the left. 
The solution f~ (x) is analytic throughout the finite plane except at the points 
congruent on the right {at a distance of n units or more) to the singularities 
of the functions a t (x), a 2 (x), . . . ., a n (x). Moreover, the singularities of each 
solution in the finite plane are poles provided that a 1 (x) , a^x) ,...., a n (x) 
are rational functions, and the complete set of numbers each of which is the 
order of one of these poles is a bounded set. 

2) There exists a right D-re-gion RD [left D-region LD] in which the 
solution f? (x) [fa(x)] is representable by the series (21) [(35)], this series 
being absolutely and uniformly convergent throughout any closed region 8 
lying entirely in RD [LD]. Equation (2) itself then suffices to determine the 
solution in that part of the plane not in the region RD [LD], 

3) If x approaches infinity in any way so as to remain always without 
any sector (however small) determined by two rays from zero to infinity and 
including between them the negative [positive] part of the axis of reals, then 

h(x) [l„(x)] being the function denoted by that symbol in the first equation 
in § 1. Also, 

li mv -(^M_l) = riimi;' n (fM-l)=ol, x = u + vV=l. 
•=±- I hi®) J L.=±- \l„(x) J J' ^ 

Remark. It is important for later use to point out the fact that, in so far 
as any statement in the conclusion of the preliminary theorem refers to a region 
RD [LD] or any part of it, it is independent of the nature of the coefficients 
^i(x), %(x), . . . ., <P„(#) of (9), and hence of the coefficients a x {x), a 2 (x), . . . ., 
a n (x) of (2), in the part of the plane exterior to RD [LD]. More precisely, 

*/ the functions ^(x), %(x), , V n (x) are analytic in some right [left] 

D-region and throughout such a region are represented asymptotically by the 
expansions in (10) or throughout such a region merely satisfy an inequality 
of the form 

!»»(*) I <]^H3» 

25 



198 Cabmichael: On the Solutions of Linear 

where p, is a constant, then there exists a solution ff (x) [f~ (x) ] and a right 
[left] B-region RD [LD] such that ff(x) [/;"(#)] has in RD [LD] the same 
properties as the solution ft(x) [f^(x)] given in the theorem. Likewise the 
inequalities derived above for such solutions in such D-regions remain valid 
under the less restrictive hypotheses. 

It should be observed, at this point, that the two solutions obtained above 
were gotten by the aid of certain functions h(x), l 2 (x), . . . . , l n (x), these being 
obtained from the formal expansions (5) by breaking them off at the (m + l)-th 
term. It follows therefore that the functions ff(x) and fn(%) are now to be 
treated as functions of m as well as of x ; or, more exactly, we are to look upon 
them as possibly functions of m. Ultimately we shall show that they are 
indeed independent of m (at least when m is sufficiently large), this proof 
being readily made after we have obtained a fundamental set of solutions, each 
solution being a possible function of m. As soon as this independence of m 
is proved, we shall be able to obtain readily the asymptotic character of ft(x) 
and /~ (x) . It is clear that the limits in the above theorem are sufficient for this. 

§ 3. The First System of Intermediate Solutions. 

For the difference equation (2) we shall now obtain four sets of n partic- 
ular solutions each, the solutions in each set being independent with respect 
to periodic multipliers of period 1, and therefore constituting a fundamental 
system of solutions. 

We shall first start from the particular solution f?(x) obtained in the 
preceding section. By means of the substitution 

f(x)=ft(x)g(x), (36) 

equation (2) goes over into the form 

g(x+n) +ka 1 (x)ft(x+n—l)g(x+n—l) + + ka n (x)f+(x)g(x) =0, (37) 

where & = l//f (x-\-n). If we put 

G(x) -=Ag(x), 

the preceding equation becomes 

G(x+n— 1) + B x {x) G{x+n—2) + . . . . +B n _ 1 (x) G(x) = 0, (38) 

where 

B^x) ^1 + ka^x) ff {% + n— 1), 

Bz(x) =l + ka 1 (x) ft (%+n — 1) +ka 2 (x)f? (x+n — 2), 

i 

B n _ 1 (x) = l + ka 1 (x)ft(x+n-l) + ....+ka n _ 1 (x)f+(x+l). 
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From the way in which equation (37) was found, it is clear that its formal 
power series solutions g t {x) may be determined from the relations 

gi{x) = a\ X 4l + Cn^ + £*-*+....y ' = 1 ' 2 ' ' •••'"• (39) 

by formal division in the second member for each value of i. Furthermore, 
the formal power series solutions G t (x) of (38) may be obtained from these 
by reckoning out formally the differences Agi(x), i = 1, 2, . . . ., n. Thus 
we have 

G i (x) = (^Yx^-^(^ 1 —l\(l + ^+^+....\ i = l,2,....,n—l. (40) 

Let us denote by L t {x), i = l, 2, . . . ., n — 1, the functions 

Li{x) = (^V-^^-l) (l + ^+ .... +^), » = 1, 2, . . • ., n-1, (41) 

where m has the same value as in the preceding section. 

An examination of equation (38) with reference to the first remark 
following the theorem of the preceding section is sufficient to bring out the 
fact that there exists a right D-region in which the coefficients B x (x), .... 
B n _ 1 (x) satisfy the conditions requisite for the application of that part of the 
theorem which relates to the right D-regions. In order to set up in detail the 
inequalities by which this may be proved, it is necessary to separate equation 
(38) into two members after the manner in which equation (2) was separated 
into the two members in (9), and then to employ the inequalities by which we 
have found fi(x) to be restricted. There is nothing involved in this more 
than the straightforward algebra of inequalities, and one can even see before- 
hand what must be the result of the computation, so that it is unnecessary 
to perform it. 

By means of the theorem of the preceding section we now have the follow- 
ing result: There exists a right D-region BD and a solution G?(x) of equation 
(38) such that G? (x) is analytic throughout BD, and is representable in such 
a region by a series of type (21) which is absolutely and uniformly convergent 
throughout any closed region S lying entirely in BD. If x approaches infinity 
in any way so as to remain always without any sector (however small) deter- 
mined by two rays from to oo and including between them the negative part 
of the axis of reals, then 

lim x m (^LJ-^-l) = 0. 
Also, 

M^f^M-lUo. 



(42) 
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These limits may be replaced by the somewhat stronger inequalities by which 
they were derived, namely, those corresponding to (31'). Then for appropri- 
ately chosen constants K and 8 we have 

\Qt{x)—L 1 {x) | < \L 1 (x)\-K\x\- m -\ «>0, 
\Qi{x)-L 1 {x)\<\L 1 {x)\-K\v\- m -\ \v\>8, 

where x=u + V\/ — 1 and lies in the region RD. 

Now, equation (37) has particular solutions g?(x) such that 

&gt{x) = Gf{x). (43) 

There are two solutions of (43) from each of which we are led to interesting 
particular solutions of (37) and hence of (2). Of these we shall treat one in 
detail in this section and the other in a later section. 
Equation (43) has the obvious formal solution 

g+( x )=Gt(x-l)+Gt(x-2)+Gt(x-3) + (44) 

In order to treat conveniently the matter of convergence of this series, it will 
be necessary to put on the further restriction that |a 1 |<|a 2 |. (Compare 
equation (20).) Then, from the second inequality in (42), it follows readily 
that the series in equation (44) converges for every value of x such that x — 1' 
x — 2, x — 3, .... all lie in RD; that is, it converges for every x such that 
| v J >Ri , where R x is the radius of the semicircle by means of which the region 
RD is denned. Moreover, if 8 is any closed region throughout which | v | >ii x , 
then the series in (44) is clearly uniformly convergent in 8. But each term 
of this series is analytic in 8. Hence, gf (x) is analytic in 8 ; it is therefore 
analytic at every point of the finite plane for which |v|>i?i. 

We may also obtain an upper bound to the absolute value of gf (x) when 
I v I =^i • We observe that R x should be chosen greater than 2 K and greater 
than 8 so as to make valid all our inequalities. Define a function 8(x) by 
means of the series 

8(x) =L 1 (x—l) +L 1 (x—2) +L 1 (x—3) + (45) 

It is easy to see that this series is absolutely and uniformly convergent when 
x is confined to a closed region S throughout which | v \ >R t . Then, from (42) , 
(44) and (45) we have 

CO 

\g+(x)-S(x) | < 2 \L t {x-i) | • K\v\- m ~\ |v|>B a . 
Now, positive constants N x and 2V' 2 exist such that 



#i 



a t 



\Li(x)\ 
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gt(x)-S(x) | <N 2 K\L 1 (x) | |^pi 5 



From (45) we see that \8(x) \/\L 1 (x) | is bounded away from zero.* From 
this fact and the last relation it follows readily that a positive constant N s 
exists such that 



S~(x)~ 



< 



N, 



»+i ■ 



(46) 



Referring to the definitions of L x {x), l x {x) and h(x), we see that 

£ > w =$SiH£j + ft) '"~&8+£i+ •-)• <46a) 



where y m+1 , . . . . are constants, the expansion being valid for sufficiently large 
values of \x\. From this relation and (45) it follows that 



S(x) 



l 2 (x) 



a a 



(x — i) f 



Ym+l 



(45') 



k(x) ' ,~W v ~ "' \(x-i) m+1 ••• 

On transposing the term Z 2 (rr)/^(a;) and dividing the resulting equation 
through by this term, we have a relation from which it is easy to show that 
a positive constant N t exists such that 



8(x)l 1 (x) 
l 2 (x) 



<X 



<x 2 



N, 



\x — % 



m+l 



N, °° 



vl m+1 i= l 



.), 



In showing this, it is necessary to observe that 

h(x) w 

the quantity in parenthesis being a descending power series in x. Hence, 
a positive constant N 5 exists such that 



8(x)l 1 (x) 
l 2 (x) 



< 



N 6 



|m+l' 



(47) 



Multiply relations (46) and (47) together member by member and combine 
the resulting inequality with (46) and (47) ; thus we have 



gH x )h(«>) 
l z (x) 

where N e is a positive constant. 



< 



N, 



6 



I m+1 > 



(48) 



* Or this may more readily be proved by aid of relation (47) and the fact that Li(x)h(a;) /h(&) 
approaches (02/01) — 1 as a approaches infinity, the latter fact being readily demonstrated by aid of equa- 
tion (46a). 
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Now, a solution f£(x) of equation (2) is given by the relation 

It (*)= ft {<*)&{*), 

as one sees from (36). This solution is analytic for \v\ sufficiently large. 
From (31') we have a relation of the form 






< 



N 7 



I m+l > 



where N 7 is a positive constant. From this relation and (48) it follows readily 
that a constant N exists such that 



l 2 (x) 



< 



N 



I m+l > 



\v\>R, 



(49) 



where R is a positive constant. This is one of the fundamental inequalities 
for Jt (%). 

A similar one may be found corresponding to the first inequality in (42). 
It may be written 






< 



,w 



\x 



m+l '• 



«>0, \v\>B, 



(50) 



where (i is a positive constant. In view of (49) it is obviously sufficient to 
prove this for the case when x is further limited to lie in a sector made by two 
rays from zero to infinity and lying entirely in the right half-plane. Further- 
more, it is clear that \x\, and hence u, may be taken as large as one pleases. 

We start from equation (44). Choose x so that u is large and denote by 
k the greatest integer such that 2k <u. Write (44) in the form 



g+(x) = ZG+(x-i)+ 2 Gf{x-i). 

i=l i=k+l 



Then 



g+(x)-8{x) = S l \Gt{x-i)-L 1 {x-i)\+ 2 \Gt{x-i)-L 1 {x-i)\, 

i=l i=k+l 



so that, in view of (42), we have 
\g}(x)—S(x)\<^\L 1 (x—i) 



K 



since it is clear that the first inequality in (42) is applicable to the first 
summation and the second inequality to the second. 

Now, it is clear that a positive constant (i x exists such that 



L t (x — i) 
Then from (51) we have 



<(l! 



i = l,2, , *. 
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2K * 



\gt{x)— S(x) | < ft |L 1 (aj) iT-r ro+1 . =f 



+ ft 1 \L 1 (x) 



K 



.,v m+1 i= 7+i 



<*2 
«ll 



If | a? | increases indefinitely, it is clear that the second term in the second 
member becomes indefinitely small with respect to the first term ; to see this 
it is sufficient to observe the role of the exponential quantity | a-j/aj | - *. There- 
fore, a positive constant /u 2 exists such that 

\gt{x)S(x)\<\L 1 {x)\- ( t t \x\-~- i . 
From this point forward the march of the argument is like that in the 
preceding treatment. As in the previous case it may readily be shown that 
\S(x) l/lLxix) | is bounded away from zero; hence we see from the last ine- 
quality that a positive constant (i s exists such that 

9\ \' X ) 1 ^ ,. l-l-m-l 

By a method similar to that by which (47) was proved, it may be shown 
that a positive constant (U 4 exists such that 

In the demonstration of this result it will be necessary to separate the series 
for S(x) in (45') into two series, one of k terms and the other of the remaining 
terms. The argument is similar to that just employed; it may readily be 
supplied by the reader. 

Making use of the last two relations and continuing precisely as in the 
argument following (47) one may readily complete the proof of (50). Thus, 
we have the requisite properties of the second solution ft{%) of equation (2) 
in case \a x \ < | a 2 1 . 

In case | a x | = | a 2 1 the method by which ft(x) was found is obviously valid 
for a second solution f£{x). We now call it f£(x). It has the properties 
expressed in (49) and (50). 

Now, apply to equation (38) the result which we have just derived for 
equation (2). "We see that (38) has a second solution G£(x) such that 

\G}(x)-L 2 (x) | < \L 2 (x) l-j^lal— \ u = 0; 
\Gt{oo)-L 2 {x) | < \L 2 (x)\-K 1 \v\- m -\ 

where K x is a positive constant, these relations being valid when | v \ is suffi- 
ciently large. If we start from this solution and proceed by a method similar 
to that employed in the discussion associated with G?(x), it is clear that we 
shall be led to a third solution f£{x) of equation (2), that this solution will be 
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analytic for | v I sufficiently large and that it will verify the following relations : 



tlM. 

l s (x) 

/7(g) 

l 8 (x) 



< 



ft' 



\x 



m+1 ' 



t*>0; 



< 



I m+l ' 



(52) 



when | v | is sufficiently large, fi' being a constant. 

By continuing this process of interaction between equations (2) and (38), 
we shall be led finally to a set of n solutions, 

ft{x)=n{x), Jt(x), .... ,/+(*), (53) 

of equation (2), each solution being analytic when \v\ is sufficiently large. 
From the inequalities such as (52), by which these functions are bounded, we 
have the following relations : 

r ft(v) 



lim x m 



lim v m 



KM 

k{x) 



-lj=0, »=1, 2, . . 
-l|=0, i = l,2,.. 



n, 



(54) 



the first relation being valid for x approaching infinity in any way so as to 
satisfy the two conditions that \v\ is sufficiently large and that x remains 
without a sector formed by two rays from zero to infinity and including 
between them the negative part of the axis of reals. 

That these solutions are linearly independent with respect to periodic 
multipliers of period 1, and hence constitute a fundamental set of solutions of 
(2), is easily shown. It is sufficient to prove that the determinant 

Jj(x) /+(*) ..../7(*) 

/i + (a + l) / 2 + (^+l) .... / + (£ + !) 



D(x)^ 



f+( x + n -l)f}(x+n—l) . . . .fi{x + n-l) 
is not identically equal to zero. Divide the &-th column of this determinant by 
af#"*, k=l, 2, . . . ., n. It is sufficient to prove that the resulting determinant 
is not identically zero. But in view of (54) it may be seen that this determi- 
nant approaches the value D ( oo ) as | v \ approaches infinity, D ( oo ) having the 
same meaning as in the first part of § 1. This is different from zero, since 
a 1? . . . ., a„ are all different. Therefore, the set (53) is a fundamental system 
of solutions of equation (2). 

We have seen that the function ft{x) was possibly dependent on m. The 
other functions in (53) were obtained by aid of ft(x), and therefore they are 



1=0. 
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possibly functions of m. We shall now show that they are indeed independent 
of m, at least if m is sufficiently large. 

More generally, let f k (x) be any solution of equation (2) such that 

lira £c m fAl^) 1) = (55) 

when x, =u-\-vV — 1, approaches infinity in such a way that u approaches 
+ oo , and | v | becomes sufficiently large, say | v | > V. 

Since f h (x) is a solution of equation (2) it may be written in the form 

f k (x) = Pl (x)ft(x) +p,{x)ft{x) + . . . . +p % (x)fi(<c), (56) 

where Pt{x), p 2 (x) , . . . . , p n (x) are suitably determined periodic functions of 
x of period 1. Then, from (55) and the first equation in (54), we see that 

ta^W «§- + *(„«£ + + p. ( ,)M| 

Let x , =^ +t; V — 1, be any value of x such that |« |> I 7 , let t be a vari- 
able running over the set 0, 1, 2, 3, ... ., and let s be a non- negative integer. 
Then on replacing x by x -\-t-\-s in the preceding limit and reducing, we have 
the relations 

lim *"|r H M£ 1 **-*> • ulAl + eu) + ■•■■ + r kn < n t^« • «J,(l + } =0, 

where % ( =a,/a fc , e u approaches zero as t approaches infinity, and 

r kl =Pi(x )ui\, lj=k, r kk =p k (x )—l. 

It is to be observed that no two of the quantities u kl are equal and that each of 
them is different from 0. 

Now suppose that m is not less than the real part of any of the differ- 
ences (i t — fi k , (i 2 — (i k , . . . . , (i n — n k . Then from the preceding limit we have the 
following : 

lim\r kl uy^ • <i(l + + • • • • +r tm «M H • «U(l+0 I =0, (57) 

where the real part of v, is zero or positive for every value of I. This rela- 
tion is valid for any non-negative integral value of s. 

Consider the set of relations obtained from (57) by giving to s the values 
s = 0, 1, 2,...., n — 1. By multiplying these in order by the cofactors (in 
order) of the elements in the l-th column of the determinant 

1 + % l+£ 2 o l+e„ 

%i(l+<?n) u k2 (l + e 21 ) u kn (l + e nl ) 



A= 



26 



««T 1 (l + ei.^-i) 
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and adding the resulting relations, we have 



Since 



lira r kl uUt v > A = 0. 

t=oa 



1 1 1 



lim A = 



we have 



U kl %2 • • • • U kn 
u kl U k2 • • • • U kn 



U kl U k2 • 



• u kn 



lim r u uiit v ' = 0. 



*0, 



This relation can be true only if \u kl \ <1 or r u = 0, since the real part of 
v t is not less than 0. But \u u \ = \aL,/a k \> 1 unless l<k. Therefore r kl = if 
l>k. Hence, p k (x )=l and p t (x )—0 if l>k. Therefore, we see from (56) 
that 

f k (x)=p 1 ( X )f + (x) + ....+p k _ 1 (x)tt_ 1 (x)+Jt(x), \v\>V, (58) 

where P\{x), . . . ., p k ^ 1 (x) are suitably determined periodic functions. 

From the particular case when k = l we see that a solution f 1 (x) of equa- 
tion (2) which satisfies the relation (55) for k = 1 coincides with /f (x), at 
least if m is greater than the real part of each of the quantities /x 2 — jUj , ju 3 — jUj , 
. . . .,(i n — (ij_. From this it follows readily that each of the solutions of (2) in 
the set (53) is independent of m, provided that m is greater than the real part 
of every one of the differences (i t — (i jt j=l, 2, . . . ., n, i=j-\-l J . . . .,n. In 
order to see this, it is only necessary to observe the manner in which these 
solutions were obtained. 

It should be observed that the only use made of the fact that the real part 
of v l is not less than zero was to insure that u\iV l does not approach zero, for 
t becoming infinite, when I > k. This condition is obviously satisfied without 
regard to the value of v t , provided that the roots a x , a 2 ,••••, a ffi of the char- 
acteristic algebraic equation are different in absolute value, for then \u kl \ > 1 
when l>k. Hence we may take m—0 for the case when the a's are all differ- 
ent in absolute value. 

We shall now determine largely the asymptotic character of the solutions 
(53) under the hypothesis that m is greater than the real part of each of the 
differences ^ — ft, j = l, 2, . . . ., n, i—j + 1, . . . ., n, or without restriction on m 
when the roots of the characteristic algebraic equation are different in absolute 
value. 
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In this connection it is convenient to distinguish (cf. Birkhoff, loc. cit., 
p. 248) between two kinds of asymptotic representation of a function by a 
series of the form 

Let g(x) be the given function; if for each m the difference 

d(x) =g{x)x-°*a-*x-> t -(c + -^ + .... + ^) 

becomes uniformly small of order m ; that is, if lim x m d(x) =0 for x approach- 
ing infinity in a certain region, we say that g(x) is asymptotically represented 
by S(x) in that region, with respect to x; if, on the other hand, we have 
lim v m d(x) =0, we say that g(x) is asymptotically represented by 8(x) in that 

v=± m 

region, with respect to v. 

Making use of relation (54) for increasingly large values of m and remem- 
bering that ft{x) is independent of m, we see readily that fti x ) is asymp- 
totic to 



a$aP 1+-JL + -S + .... 

\ X X 1 

with respect to v and also with respect to x, provided that x approaches infinity 
so as to satisfy the conditions that \v\ is sufficiently large and that x remains 
without a sector formed by two rays from zero to infinity and including between 
them the negative part of the axis of reals. 

The fundamental properties of the functions fi(x), . . . ., fn(x) are now 
determined throughout a region in which | v \ is sufficiently great, say | v | >F, 
and x is to the right of some line parallel to the axis of imaginaries. Their 
properties in the remainder of that part of the plane for which | v | > V may be 
determined, as in the case of ft{x) in § 2, by means of equation (2) itself. 

The principal results indicated in this section may be put together in the 
form of the following theorem : 

First System of Intermediate Solutions.* If the characteristic algebraic 
equation associated with the difference equation (2) has its roots different from 
each other and from zero, then the difference equation has the fundamental sys- 
tem of solutions 

♦The existence of these solutions was first pointed out in my memoir, loo. cit., p. 119. Their prop- 
erties were first developed in Birkhoff's memoir, already referred to. 
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fi{^),li{x),....Ji{x), (53 bis) 

determined above, and these solutions have the following properties: 

1) A constant V exists such that each function in the set is analytic when 
\v\>V, where v is defined by writing x—u+vV — 1, u and v being real; 

2) The asymptotic relations 

ft{*)~w(l+^ + % + ....)> » = 1, 2, . . . . , n, 

\ x x I 

exist and are valid with respect to x in that part * of any right half-plane for 
which \v\>V and with respect to v in the entire plane. 

(By a right [left] half-plane is meant that part of the plane which is to 
the right [left] of a line parallel to the axis of imaginaries.) 

Let us retain for m the value if the roots of the characteristic equation 
are different in absolute value and in the contrary case a value not less than 
the real part of any of the differences fi t — /u,-, j = l, 2, . . . ., n, i—j-\-l, . . . ., n. 
Denote by li(x), i=l, 2, . . . . , n, the functions so represented in the first part 
of § 1. Then a solution f(x) of equation (2), 

f{x)=p 1 {x)li{x)+p i {x)]t{x) + ....+p n {x)]Z{x), (59) 

where p^x), p 2 (x), . . . ., p tt (x) are periodic functions of x of period 1, has an 
important set of properties obtainable from the first relation in (54), namely, 



t=x 



limtm ( f(xo+t)-pMfn%o+t) teu-o i«i>f 



t=<a 



lirv, ^.f /foo+O — P n(^o)fn(Xo+t)-.....~P2 (Xo)f}( ^0+t) _ / _ y { n 



t = co 



(60) 



The proof of this is so far similar in character to the argumentation associated 
with equations (55) to (58) that it may readily be supplied by the reader. 
Furthermore, any solution of (2) having the properties (60) is easily seen to 
be identical with the solution f(x) defined by (59). Consequently, the relations 
(60) may be looked upon as a sort of set of "initial conditions" at infinity, 
defining uniquely a solution of the given equation (2). 

In general these relations are complicated. It may be observed, however, 
that the single " initial condition " 

♦The asymptotic character of the first solution ft(®) is maintained with respect to x in any unre- 
stricted right half-plane. 
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l rAw>-"^}=° 



is sufficient to define uniquely a solution f(x) =p(x)ft{x). In this case alone 
is the solution defined in a simple way by the initial conditions. 

§ 4. Second System of Intermediate Solutions. 

The fundamental set of solutions of the preceding section was obtained by 
means of a class of transformations depending on the solution ft(x) of § 2. 
If one starts from the solution fa(x) of § 2 and proceeds by a method in every 
respect similar to that employed in § 3 he will obviously be led to a new funda- 
mental system of solutions related to those in § 3 as f„(x) is related to fi(x). 
It is unnecessary to give the argumentation. We merely state the result : 

Second System op Intermediate Solutions. If the characteristic alge- 
braic equation associated with the difference equation (2) has its roots different 
from each other and from zero, then the difference equation has a fundamental 
system of solutions, 

7r(«),7r(«), — ,/»(«), 

possessing the following properties: 

1) A constant V exists such that each function in the set is analytic if 
\v\>V, where v is defined by writing x=u-\-vV — 1, u and v being real; 

2) The asymptotic relations 

/r(z)~aWl+ ^ + ^ + . . . .), i = l, 2, ..... n, 

\ {JO 30 / 

exist and are valid with respect to x in that part* of any left half -plane for 
which | v | > V and with respect to v in the entire plane. 

The second property implies a set of relations analogous to those in (54). 
In fact, it is by means of such relations that this part of the theorem is demon- 
strated. From these follows a set of " initial conditions " analogous to those 
in (60), the limits in the present case being taken for t approaching — oo . The 
reader may readily supply them. 

§ 5. The First System of Principal Solutions. 

By means of the particular solution ff(x) of equation (2), obtained in § 2, 
we proceeded in § 3 to build up a fundamental set of solutions. In doing this 

♦The asymptotic character of fa(a>) is maintained with respect to x in any unrestricted left half- 
plane. 
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we made a transformation which introduced a new function gt{x) satisfying 
the relation 

Ag+ (x) =Gf (x) , (43 bis) 

Gf (x) being a particular solution of equation (38), this latter equation being 
obtained from equation (2) by means of the transformation (36). For g?(x) 
we took the function gf(x), defined in equation (44), this function satisfying 
equation (43). We stated that another solution of (43) would also lead to 
interesting solutions of (2). These latter solutions we shall now find. 

Throughout this discussion we shall assume that the integer m has been so 
chosen that ff(x) is independent of m, this being possible, as was shown in § 3. 
Then we have not merely the limits affecting ft(x), which were stated in the 
preliminary theorem of § 2, but also the further property expressed in the 
asymptotic relation 



\ X x I 



valid with respect to x in the entire right half -plane and with respect to v in 
the entire plane. Similarly, if we make use of this extended property of ft (x) 
then, instead of the inequality (42), we have the asymptotic relation 

valid with respect to x in the entire right half -plane and with respect to v in 
the entire plane. Furthermore, the function gf(x), defined by equation (44), 
is asymptotically represented by 

as one sees readily from its properties as developed in § 3 ; the range being the 
same as before, except that | v | must now be restricted to be not less than a 
certain fixed quantity. 

A solution gf{x) of (43) in the form of a contour integral, the path of 
integration going to infinity along two parallel lines, will be found to serve our 
purpose. We form the function 

gi(v)=SLp(x— z)p(x, z)Gt{z)dz, (61) 

where 

p (x— z) = \ 1 — e **V-K.-,) | -i ; 

and where the path of integration L is yet to be chosen and the function p (x, z) 
is to be determined subject to the condition that it is analytic with respect to x 
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and also with respect to z throughout the finite planes of these variables and 

that 

p(x+l, z)=p(x, z). 

Let RD be a right D-region throughout which ft (x) and G} (x) are both 
analytic and verify the relations (31') and (42), respectively, for every value 
of m. Since G?(x) is analytic throughout RD, it is clear that all the infinities 
of the integrand in (61) are at the points z=x±r, where r is zero or a posi- 
tive integer. Let AB be the straight line containing all these points. Let the 
path of integration L or CKRH lie entirely within RD and be formed in the 
following manner : It consists of three parts : a part RH lying entirely above 
AB and extending to infinity in a negative direction RH parallel to the axis of 
reals ; a part KG lying below AB and extending to infinity in a negative direc- 
tion KG parallel to the axis of reals; a finite part KR which crosses AB 
once between x — 1 and x and at no other point. Moreover, the path CKRH is 
such that the part of the finite plane not in RD lies entirely to the left of this 
path. 

In order that the integral along such a path shall exist, it is necessary 
that the function p (x, z) shall have certain requisite properties. The value 
p{x, s)=l would satisfy this requirement, provided that |a 1 |<|a 2 |; and it 
would indeed lead to solutions of our equation (2) (compare my memoir, loc. 
cit., p. 119). Simpler solutions may, however, be obtained by a different choice 
(compare Birkhoff, loc. cit., p. 264). 

For the present we assume that |a x | < |a 2 | and choose for p(x, z) the 
function 

p(x, z)=e 2rxy - l(x - e \ 

where X is the least integer not less than the real part of 

2 -^= f (loga 2 -lo g « 1 ). 

Then for gt{x) we have the value 

g+(x) =f L p(x—z)e 2 ^- lix -°>G+(z)dz. (62) 

It is evident that this function is analytic throughout RD. Furthermore, the 
value of gt (x), for a particular x, is unaffected by any deformation of the path 
of integration, provided that during the deformation it has always the proper- 
ties stated in its definition. 

We show that gt{x), as defined in (62), is indeed a solution of equation 
(43). We have 

gt (x+ 1 ) = f v p (x-z) e^V-iAo,-)^ {z) dZ} 
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where L' is a path CKTSRH which differs from CKRH only in crossing AB 
between x and x-\-l instead of between x — 1 and x. Then 

J h< J L = J KTSRK • 

By Cauchy's theorem on residues the last integral is equal to 2tiV — 1 times the 
residue of the integral at the point s=x; that is, it is equal to Gt(x). Hence, 

gt{x + l)-gt{x)=Gt{x); 
or, g x (x) satisfies equation (43). It is this particular solution of (43) which 
we shall denote by g?(x) in the remaining discussion. 

We shall now determine a bound to the magnitude of \yt {%) \ for certain 
ranges of variation for x. For this purpose it is convenient to break up the 
path of integration into two parts L x and L 2 , as follows : The contour L x is the 
fixed contour oo A x B x oc lying entirely in the region RD and containing within 
it that part of the finite plane which is exterior to RD. Moreover, it consists 
of three parts, each of which is a straight line, a finite part A X B 1 perpendicular 
to the axis of reals, and two infinite parts A x <x> and B x <x> parallel to the axis of 
reals. The contour L 2 lies entirely without the region inclosed by L x . If x lies 
above B x <x> , L 2 consists of a loop circuit to infinity lying above B x <x> and includ- 
ing within it the points x — 1, x — 2, x—3, . . . ., but not the points x, x+1, x + 2, 

Moreover, in the neighborhood of infinity the bounding lines of L 2 are 

parallel to the axis of reals. If x lies below A x cc , L 2 consists of a loop circuit 
to co of the same nature as when x lies above B x <x> . If x lies between B x x> and 
Ax<x> , then L 2 consists of a finite loop containing within it the points x — 1, x — 2, 
. . . ., x — I, but no other point congruent to x, x — I being the leftmost point 
congruent to x and not on or within the contour L x . For the integrand in (62) 
it is clear that the integral along L is equal to the sum of two integrals, one 
along L x and the other along L 2 , provided that x is not congruent to a point on 
the contour L x . In symbols 

where I L denotes the integral in (62) and I Ll and l Lz similar integrals taken 
along the paths L x and L 2 , the integrand being the same in all three cases. 

Let us consider first the integral I Li It denotes a periodic function of x 
of period 1 since the contour is fixed and the integrand is a periodic function 
of x of period 1. It may be written in either of the forms : 

£ i J Ll 1 e 2irV-10c-z) ' 

t _A(x-i)V-i.| ? "i vgj as 

Ll Jl l 0-27TV-K3-3) ^ 
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If x, =u + vV — 1, approaches infinity so that v becomes positively infinite, the 
integral in the first one of these equations is bounded. If v becomes negatively 
infinite, the integral in the second one of these equations is bounded. If x 
approaches infinity while v remains bounded, then I Li is itself bounded. Hence, 

'e 2T ^~ lx q(x), 

e^ x - l ^~ lx q(x), 

according as x lies above -Bjoo , between B x cc and A^ , or below A-^cc , the 
function q(x) in each case being periodic in x of period 1 and bounded. This 
result is strictly valid only when * is not congruent to a point of the contour 
L x . For all such x as this, it is obvious that a single modified contour might 
be used and that a similar form for I L would be obtained, the periodic functions 
q(x) being different in this case, but still being bounded. 

Let us next consider the integral I Ll . We separate the treatment into 
two cases. 

In the first place, let us suppose that x lies above 5 x co or below Aioo . In 
this case the integral l Li may be evaluated as a sum of residues multiplied by 
2tiV — 1 ; thus, we have 

I Lt =Or(x-l)+Gf(x-2) +Gt(x-3) + 

Hence, in this case, 

a function whose asymptotic character is known. 

In the second place, let us suppose that x lies between 5 x oo and Ai<x> and 
to the right of A X B X . Then, evaluating as before by means of residues, we have 

I L = kGt(x-k). 

The asymptotic character of I Ls with respect to x in the right half -plane can 
now be determined in precisely the same way as that by which the correspond- 
ing bounds to the increase of \gf (x) | were determined in § 3. The only modi- 
fication necessary arises from the fact that we now have a finite series instead 
of the infinite series in (44), so that when this series is separated into two parts, 
as in the earlier discussion, the second part as well as the first consists of only 
a finite number of terms. But it was precisely this second part which did not 
affect the asymptotic form, so that the modification in the argument will not 
affect the conclusion. 
27 



214 Cahmiohael: On the Solutions of Linear 

Taking these two cases together, then, and remembering the asymptotic 
form of gt(oo), we find that 

I ~^V-M1+ ^ af^(l + c 21 a;- 1 +....) 
L ° W V -l----^ a >af*(l + c ll ar 1 +....) 

throughout the right half -plane. 

Now let us define a second solution ft{oo) of equation (2) by means of 
the relation 

ft{v)=9t{x)ft{x) = {l Ll +l Ll )ft{x). 

From the known form of l Ll and the known asymptotic character of I L2 and 
/J (x) we see that 



(63) 



fi{x)~e i ^- lx q{x){a\% h (l+ ^ + . . . \\ + a * x »*(l+ c Jk + 

ft{x)~q (»)[«;af»(l+ — + ■■■ .)]+a&*(l+ °f + . . . .), 

ft(x)^e 2 ^ x - 1)y/ - u q(x)Lfx^ (l+ °* + . . . . Yl + a2 Wl + ^ + . . . . V 

according as x lies above j^co , between -Sjoo and A x <x> or below J. a ao , the 
asymptotic representation being valid with respect to x in any right half -plane 
and with respect to v in the entire plane. 

From this result we may determine, in simpler form, the asymptotic char- 
acter of ft(x) with respect to a; in a right half -plane. Suppose first that x 
lies above i^cc . The dominating term in the above representation of ft (x) 
depends upon the dominating term of the set 

g(27rXV-i+loga,).i; gOoga^s^ 

If we divide both of these quantities by the last, the exponents are 

log Oi— log <x 2 



. r log Oi— log <x 2 T 

2W=i{*+ 27tV==I -)«, 



0. 



If the real part of the first exponent is negative the latter of the two quantities 
predominates. From the definition of X and the fact that we now have | a x | < | a 2 1 
it follows readily that the real part of this first exponent is indeed negative 
if u (where * = M+-yV- 1) is positive and also if u is merely bounded below j 
provided that the real part of (log <x 2 — log a 1 )/(27tV— 1) is not the integer A 
and v is sufficiently large. Therefore, 

ft (x) ~<xf a* (l+^+ • • . • ) (64 ) 
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for x above B x <x> and to the right of the axis of imaginaries ; and, unless the 
real part of (log a 2 — log a,i)/ (2n V — 1) is an integer, the asymptotic represen- 
tation is valid for x above B 1 co and to the right of any line parallel to the axis 
of imaginaries. 

Suppose next that x lies below A x <x> . The dominating term in the asymp- 
totic representation (63) of ft{x) now depends upon the dominating term of 
the set 

The second of these predominates, provided that the real part of 

2w~iU-i+ l0gai ~i^ a2 U 

I 2tiV— 1 J 

is negative ; and this real part is in fact negative if u is positive or if u is 
merely bounded below and — v is sufficiently large. Hence, the asymptotic 
representation (64) is valid for x below J-xoo and to the right of any line 
parallel to the axis of imaginaries. 

Finally, suppose that x lies between Ai<x> and B 1 oo . Then the dominating 
term in the asymptotic representation (63) of ft{%) is clearly the last term, 
since now | oci | <C _j oc^ | , so that the asymptotic representation (64) is valid for x 
between -4 x co and ^oo . 

Combining these results, we see that the asymptotic representation (64) 
is certainly valid with respect to x for x approaching infinity in any way so as 
to remain to the right of the axis of imaginaries when in the upper half -plane 
and to the right of any line parallel to the axis of imaginaries when in the lower 
half -plane. Moreover, if the real part of (log a 2 — log a x )/ (2 n V — 1) is not an 
integer, then the asymptotic representation is valid with respect to x in any 
right half -plane. 

The above result was obtained on the supposition that ■ | tn t | < | ct 2 J . If 
|a 1 | = |a 2 |, we form a second solution f£{£) in the same way as ft(x) was 
formed in § 2 (see the preliminary theorem), this being obviously possible 
since a x and a 2 then play exactly parallel roles. This solution f£(x) has the 
asymptotic representation (64) with respect to x in any right half -plane. 

Thus, in any event, we have two solutions of equation (2). We proceed 
further as in § 3. Equation (38) is now known to have a second solution and 
its asymptotic character is determined in the same way as that of f£(x). By 
means of this second solution of equation (38) we are led to a third solution 
of equation (2). Then by a further reaction between equations (2) and (38) 
we find another solution of (2), and so on until n solutions of (2), are found. 
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These are analytic in a suitable right half-plane and are represented asymp- 
totically by the formal power series solutions of (2). By means of their 
asymptotic form it may be shown (as in the case of the solutions obtained 
in § 3) that they are linearly independent with respect to periodic multipliers 
of period 1, and hence form a fundamental system of solutions of equation (2). 
Since they are analytic in a suitable right half-plane, the possible position of 
their singularities may be determined as in the case of ft (*) in § 2. Combining 
the various results which, we thus have, we obtain the following theorem : 

First System of Principal Solutions. If the characteristic algebraic 
equation associated with the difference equation (2) has its roots different from 
each other and from zero, then this difference equation has the fundamental 
system of solutions 

ft(x),ft(a>),....,f:(x), 

determined above and possessing the following properties : 

1) Each function in the system is analytic throughout the finite plane, 
except at the singularities of the functions (34) and points congruent to them 
on the left. Moreover, the singularities, in the finite plane, of each function in 
the solution are poles, provided that %(«), a 2 (x), . . . ., a n (x) are rational func- 
tions and the complete set of numbers, each of which is the order of one of 
these poles, is a bounded set. 

2) The functions ft(%), *=1, 2, . . . ., n, have with respect to x the asymp- 
totic representation 



A + (s)~a«Wl+ ^ + % ■+ . . . .), i=l, 2, 

\ x x I 



n. 



this representation being valid for x approaching infinity in any way so as to 
remain to the right of the axis of imaginaries when in the upper half-plane 
and to the right of any line parallel to the axis of imaginaries when in the 
lower half-plane* Moreover, this representation is valid in any right half- 
plane provided that no one of the quantities 

log a t — log a, • -, r, ■ • . -, 

27iV — 1 

has its real part equal to an integer. 

* The upper and lower half-planes enter here unsymmetrically. Their r61es may be interchanged 
by a different choice of the quantities X entering into the discussion. See the paragraph containing 
equation (62). 
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From results previously obtained we know, furthermore, that ft{x) always 
maintains its asymptotic character with respect to x in any right half -plane and 
with respect to v in the entire plane. 

The solutions mentioned in the foregoing theorem are called principal solu- 
tions on account of the important properties which they possess. 

If f k {x) is any solution of (2) which is analytic in some right half -plane 
and is there represented asymptotically with respect to x by 

\ X XT 

then it may be shown (as in a similar case in § 3, beginning with equation (55) ) 
that* 

f k (x) =p 1 {x)ft{x) + .... +p k . 1 (x)fL 1 (x) +ft(x), 

where p Y («),... ., p k -i(x) are suitably determined periodic functions of x of 
period 1. By means of this result one again has a theory of " initial conditions " 
analogous to that associated with equation (60). 

§ 6. The Second System of Principal Solutions. 

In the preceding section a fundamental system of solutions was obtained 
by starting with the solution ft{x), obtained in § 2, and carrying out a certain 
process of transformation of equation (2) and interaction between it and the 
transformed equation. Similarly, one may start from the solution fn(%)> ob- 
tained in § 2, and by a similar process arrive at another fundamental system 
of solutions of equation (2). We state merely the result of this work in the 
form of the following theorem : 

Second System op Principal Solutions. If the characteristic algebraic 
equation associated with the difference equation (2) has its roots different from 
each other and from zero, then this difference equation has a fundamental system 
of solutions 

fT{x),f^{x), ,fn{x), 

possessing the following properties : 

1) Each function in the system is analytic throughout the finite plane 
except at points congruent on the right (at a distance of n units or more) to 
the singularities of the functions a t (x), a 3 (x) ,...., a n (x) . Moreover, the 
singularities of each solution, in the finite plane, are poles provided that a x (x) , 

* If /*(«) is replaced by the fc-th solution of the first fundamental set of associated solutions, then 
it is easy to see that we continue to have the same relation as that given in the text, this relation cer- 
tainly being valid if \v\ is sufficiently large. 
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a 2 (x), . . . ., a n (x) are rational functions and the complete set of numbers, each 
of which is the order of one of these poles, is a bounded set. 

2) The functions fT {%), i=l, 2, . . . ., n, have with respect to x the asymp- 
totic representation 

fir(aO~a?Wl+ - + -z + • • • •)' »=1, 2, . . . ., n, 

\ x x / 

this representation being valid for x approaching infinity in any way so as to re- 
main to the left of the axis of imaginaries when in the lower half -plane and to the 
left of any line parallel to the axis of imaginaries when in the upper half-plane* 
Moreover, this representation is valid in any left half-plane, provided that no 
one of the quantities 

log a*— log a, • • , i 

has its real part equal to an integer. 

Further remarks, similar to those at the close of the preceding section, 
may be made with reference to the solutions fr(%), ■ ■ • •> fn( x )- 

§ 7. Other Results. 

The methods by which the foregoing results are demonstrated are essen- 
tially distinct from any at present existent in the literature of difference equa- 
tions. There are three other results which are necessary to complete the first 
fundamentals of a general theory of equation (1). These are due to Birkhoff 
(loc. cit.), who has demonstrated them for a system of n equations of the first 
order rather than for a single equation of order n. Birkhoff's proof, with 
merely formal modifications, will apply to equation (1) and may be based on 
the theorems demonstrated above. It is therefore unnecessary to repeat the 
proof. It seems desirable, however, to have a careful statement of these in a 
form to be directly applicable to equation (1) and for the sake of completeness 
of the fundamental results assembled in this paper. Consequently, a statement 
of these theorems, without proof, is given below : 

1. If the characteristic algebraic equation associated with the difference 
equation (1) has its roots different from each other and from zero, and if the 
coefficients a a (a;), a 2 (x), . . . ., a n (x) are polynomials in x and if further we 
write 

* Compare the foot-note to the previous theorem. 
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Fr(x)= Pil (x) Ft(x)+p i2 (x) Ft{x) + . . . . + p tn (x) F+(x), t = l, 2, . . . ., n, 

where 

Ft(x)=afft{x), FT(x)=af'fT(x), i = l, 2, . . . ., n, 

then the fundamental periodic functions Pij{x) are of the form 

p a (x) = l + cff e 2 * V-i * + c£ ?) e** V-i * _|_ _|_ e trn t V-i e 27r^ V-i * f 

p. y («) = e 2,rXt > V-i * [ c $> -f- c$f e«»- V=i * _f_ + c ^-d e 2(/»-i)«- V-i *] , i zf= j f 

where 7 Hj denotes the least integer not less than the real part of 

log a t — log a, 

2*V— 1 

If the functions a x {x), a 2 {x), ...., a n (x) are restricted merely to be 
rational, then the periodic functions Pij(x) are rational functions of e 2 *^- 1 *. 

2. Let the characteristic algebraic equation associated with the difference 
equation (1) have its roots different from each other and from zero, and let 
the coefficients a^x), a 2 (x) , . . . . , a n (x) be polynomials in x. Write 



271 V ~H 2*v=i + ^j-^' 



where /l# is the least integer not less than the real part of the negative of the 
first term in parenthesis. Mark the two sets of points 

*m » e«,»-i , , e<,<-i » and <? <n — 2 rt V— 1, , <? <( i-i —2 rt V— 1, 0, 

and call them, respectively, 

Construct the convex broken line P n P a . . . . P t above which all the remaining 
points P lie, and likewise construct the convex broken line Q n Q e . . . . # t - below 
which all the remaining points Q lie. Let the acute angles which the successive 

sides of P n P ff . . . .Pi and Q n Q e Q { make with the axis of reals be $„,$„, 

and "\> n , "b e , . . . . , respectively. The critical rays between arg x—0 and arg#=2?i 
are then 

ft ft ft . | 71 , 

2 — 9*»> 2 _ ^ ff ' ' ~~ 2 *' — 2 "' 

in angular order. The asymptotic form of F^(x) changes from 

e 2rt 4 , V-i . c f X vx a- a?, M + ffi + £* + \ 



to 

e 



" \ X x 2 I 
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along the critical ray arg x—\n — ty t , and likewise from* 

e 27r(X 4 ,+« 4 ,-l)V^ia! c (/t+« 1 ,-l) e -2ir / ii,V-l a ./w a * a ./t,/-£ , £fl , j 

•' t y x • ■ ■ ■ J, 



to 

e 



iTrCK is +S i ,-l)^-\x c \li.+5i s -l) e -iirit, s ^-i x ii.x a x x ii.J-^ _j_ £«1 i \ 



along the critical ray arg#= — In+^f Between the last critical ray of the 
first set and the first of the last set the asymptotic form is given by 

a;"*a?Wl+ — + ). 

A similar theorem exists for the functions F? {%) in case equation (1) is 
such that it may be written in the form 

a' (x)F(x+n) + .... + a' n _ 1 (x)F(x + l) + F(x) = 0, 

where a' (x), a[(x), . . . ., a' n _ x (x) are polynomials 

3. Let Ff (x) ,...., F+ (x) and Ff (x) ,...., F~ (x) be two sets of single- 
valued functions which are analytic throughout the finite plane, except for 
poles, and which have the further property that 

lim Ff ( x ) x'^^x-^ = 1, lim Fr (x) x^a^x^ = 1 , 

(M being an integer, the constants a 2 , Oj ,...., a„ being different from each other 
and from zero and the range of x for the limits being as follows : In the first 
[second] case x may approach infinity in any way so as to remain to the right 
[left] of any line parallel to the axis of imaginaries when in the upper [lower] 
half-plane and to the right [left] of any line parallel to the axis of imaginaries 
when in the lower [upper] half-plane. Furthermore, let these two sets of 
functions be connected by the relations 

Fr(x)=p a (x)F?(x) + . . . .+p in {x)Ft{x), *=1, 2, . .. ., n, 

where the functions Pij(x) are periodic in x of period 1. Then the sets Ff (x), 
. . . ., F„(x) and F^~ (x), . . . ., F~(x) are fundamental systems of solutions of 
a difference equation (1) in which the coefficients a^x), . . . ., a n (x) are rational 
functions of x. 

Univebsity of Illinois. 

* Here Sij is equal to unity or zero, according as i is or is not equal to /. 



